Optical conductivity of the Holstein polaron 
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The Momentum Average approximation is used to derive a new kind of non-perturbational ana- 
lytical expression for the optical conductivity (OC) of a Holstein polaron at zero temperature. This 
provides insight into the shape of the OC, by linking it to the structure of the polaron's phonon 
cloud. Our method works in any dimension, properly enforces selection rules, can be systematically 
improved, and also generalizes to momentum-dependent couplings. Its accuracy is demonstrated by 
a comparison with the first detailed set of three-dimensional numerical OC results, obtained using 
the approximation-free diagrammatic Monte Carlo method. 
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Although the study of polarons is one of the older prob- 
lems in solid state physics [1], a full understanding of 
their properties is still missing. This is especially true for 
the excited states which influence response functions like 
the optical conductivity (OC). OC measurements have 
revealed the role of the electron-phonon (el-ph) coupling 
in many materials, e.g., cuprates [2] and manganites [3]. 
In particular, the shape of the OC curve is important, as 
it signifies large versus small polaron behavior [4]. 

The OC of polarons has been studied numerically us- 
ing exact diagonalization in ID [ ] and for small clus- 
ters in higher dimension (here, finite size effects can be 
an issue) [G]. There are also some diagrammatic Monte 
Carlo (DMC) results for the 3D Frohfich and 2D Holstein 
models [ ] . DMC gives approximation- free results in the 
thermodynamic limit, but it requires significant compu- 
tational effort; this is why there are very few DMC OC 
sets available in the literature. The conceptual problem 
associated with all numerical methods, however, is that 
they do not provide much insight for understanding the 
shape of the OC and its relation to the properties of the 
polaron. Analytical expressions are needed for this, but 
most prior work was limited to perturbational regimes [ ]. 
The one exception is work based on the dynamical mean- 
field theory (DMFT) [ ] , which however ignores current 
vertex corrections. The consequences are discussed be- 
low; here we state only that a complete understanding of 
the shape of the OC is still not achieved using it. 

It is, then, hard to overemphasize the need for an accu- 
rate analytical expression establishing a nonperturbative 
structure of the OC. In this Letter we obtain such an 
expression using a generalization of the Momentum Av- 
erage (MA) approximation. MA was developed for the 
single-particle Green's function of the Holstein polaron 
[10] and then extended to more complex models [11], in- 
cluding disorder [12]. It is non-perturbational since it 
sums all self-energy diagrams, up to exponentially small 
terms which are discarded. MA becomes exact in vari- 
ous asymptotic limits, satisfies multiple spectral weight 



sum rules, is quantitatively accurate in any dimension, 
at all energies for all parameters except in the extreme 
adiabatic limit, and can be systematically improved [10]. 

Here we show how to use MA to calculate two-particle 
Green's functions, needed in response functions like the 
OC. Besides efficient yet accurate results at any cou- 
pling, this finally provides the explanation for the physi- 
cal meaning of the shape of the OC. Moreover, this MA- 
based approach generalizes to OC calculations for models 
with momentum-dependent el-ph coupling [11]. 

We use the Holstein model [13] as a specific example 
since some numerical data is available for comparison: 




Here, cj^ and are electron and boson creation operators 
for a state of momentum k (the electron's spin is trivial 
and we suppress its index). The free electron dispersion 
Ek — —2tJ2'i^iCOs{kia) is for nearest-neighbor hopping 
on a c?-dimensional hypercubic lattice of constant a, and 
the Einstein optical phonons have energy i7. The last 
term describes the local el-ph coupling g J^i + 
written in k-space. All sums over momenta are over the 
Brillouin zone and we take the total number N of sites 
to infinity. We set ?i = 1 and a — 1 throughout. 

For the case we study here, i.e. a single polaron at 
T = 0, the optical conductivity is given by the Kubo 
formula [14]: 
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with V the volume, IV^o) the polaron ground state (GS), 
and the charge current operator j — ^et^^smqcl^c^ is 
in the Heisenberg picture. Here, e is the electron charge 
and q the component of q parallel to the electric field. 
Our main result is that the optical absorption equals: 

Re [aMA(.) (^)] = E P^'^fi'^i^)- (2) 
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Here i > is the level of MA^'^ approximation, denoting 
an increasing complexity of the variational description of 
the eigenstates [10]. However, the physical meaning is the 
same: Pii is the GS probability to have n phonons at the 
electron site, while fn\i^) are spectral functions describ- 
ing the electron's optical absorption in this n-phonon en- 
vironment. Eq. (2) shows the direct link between the 
OC and the structure of the polaron's phonon cloud. 

Further discussion is provided below. First, we derive 
Eq. (2) so that the meaning of various quantities becomes 
clear. Expanding the commutator and doing the integral 
in Eq. (1), we find cr{ui) = cr+(w) -|- (cr+(— w))*, with: 



a+iuj) = —{^o\jG{u + EoYMq), (3) 
ujV 

where G{uj) = [uj + irj — H]~^ with r/ — 0+, and Eq is the 
polaron GS energy. The usual route is to use a Lehmann 
representation, leading to the well-known formula: 

= ^ E \{i^o\Mn)fS{u; + Eo- (4) 
uj V — 



in terms of excited polaron eigenstates |V'n), En- Instead, 
we use twice the resolution of identity to rewrite: 
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xF„0(q,Q,c^ + So)(/3|cQ|Vo), (5) 

where _F'ct0(q, Q, a;) = (a|cqG(a;)cQ|/3). Since is 
the polaron GS, {|a)} and {|/3)} are phonon-only states. 
Moreover, because of invariance to translations, their mo- 
mentum is — q, respectively — Q. Eq. (5) is exact. 

Consider now these matrix elements within MA^^', 
whose variational meaning is to expand polaron eigen- 
states in the basis {cj (6])"|0)}, (V)i,j,n [10, IG]. Then, 

|a) ^ |-q,n) = ^ E» e^"*'^' (^I)"|0> > since only such 
states will have finite overlaps in Eq. (5). The sums over 
a, /? are now sums over phonons numbers n,m > 1. Note 
that n, TO = do not contribute to the regular part of 
a{uj) since | — q,0) ^ <5q,o|0), and sing • S^q — ^ 0. They 
do contribute to the Drude peak, 1)6 (w). 

The calculation of the single electron Green's func- 
tions F„TO(q, Q,w) — (— q, n|cqG'(w)cq| — Q,to) and of 

the residues (V'ojcql — q, ri) is now carried out. Details 
are provided in the supplementary material at the end 
of this article. Here, we note that because of the odd 
sing, sin Q terms, only the part of F„m(q, Q, w) propor- 
tional to ^q,Q(Sn,m has non- vanishing contribution to Eq. 
(5), explaining the single sum over n > 1 in Eq. (2). Fi- 
nally, the overlap |(V'o|Cq| — q, n)|^ is linked to the prob- 
ability P„ to have n phonons at the electron site in the 
polaron GS [17]. The P„ expressions are listed in the sup- 
plementary material. Altogether, we obtain the result of 




FIG. 1. (color online) Phonon statistics for a 3D Holstein 
polaron with (a) small to intermediate, and (b) large effective 
couplings, and Q,/t — 0.5. MA^'^^ (dashed lines, open symbols) 
and MA*-^' (solid lines, solid symbols) results are shown. The 
Lang-Firsov (LF) asymptotic result is shown for A = 1.25. 



Eq. (2), where: 
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X [Ao{q,cu + Eq - nQ) - Ao(q, -i^ + Eo - nn)], (6) 

with Ao{q,uj) — —^lmGo{q,uj) being the free electron 
spectral weight, where Go(q,w) = {uj+ir]—eq)~^ . Within 
the n-phonon sector of \4'o), the electron can carry any 
momentum q and (w) basically describes its optical 
absorption in the presence of this phonon environment. 

At the MA'^^ level, the variational basis is supple- 
mented with the states {cj(6])"6j'|0)}, (V)i,j 7^ l,n, 
which are key to describe the polaron-|-one-phonon con- 
tinuum [ ]. This enlarges the {\a)} and {|/3)} sets with 
the states |-q,n,5) = ^ e"*'*''^' (&I)"^I+5|0) for ah 
(5 7^ 0. The new matrix elements are found similarly and 
give much lengthier yet more accurate formulae for Pn^^ 
and fn\u!) [ ], but with the same physical meaning. 

We now discuss the key features of the 3D OC results. 
This is because (i) no accurate Holstein OC results were 
available in 3D, and we provide the first detailed set of 
DMC data; and (ii) DMFT is a better approximation for 
higher D. For the OC, the key approximation of DMFT is 
to ignore current vertex corrections, so we test its validity 
in 3D where it should be most accurate. 

GS phonon probabilities P„ for a 3D Holstein polaron 
at several eflFective couplings A = g'^/{2dtil) are shown in 
Fig. 1. ID results are quite similar [17]. At weak el-ph 
couplings, the polaron has few phonons in its cloud so 
Pn>i are small. At the crossover to a small polaron at 
A ^ 1, the distribution changes abruptly and dramati- 
cally. As A 00, F„ — > t^^^ Lang-Firsov 
(LF) limit, already quite accurate for A = 1.25. The 
difference between MA^^' and MA^^^^ is small. This is 
not surprising, since these probabilities are for the GS, 
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which is aheady very accurately described by MA^"^ [10]. 
The additional basis states added at the MA^^^ level are 
essential to describe excited states in the polaron+one- 
phonon continuum, starting at Q above the GS, and due 
to a phonon excited far from the polaron cloud [''']. As 
we show now, they do have a significant effect on the 
/„(a;) functions and therefore on the OC onset. 

The w-dependence of the OC is dictated by /„(a;). Eq. 
(6) shows that fn^ (w) becomes finite at uj + Eq — > 
—2dt, because the free-electron spectral weight is finite 
in [—2dt,2dt]. This implies that the onset of absorption 
is set by the n = 1 curve to be Wth = —2dt — Eq + H., and 
larger n contributions are shifted (n — 1)0. higher. As A 
increases and Eq falls further below ~2dt, this suggests 
that Wth increases monotonically. This is wrong: the OC 
onset is always expected at ojth = [ , ]• The discrep- 
ancy is easy to understand. The onset is due to absorp- 
tion into the polaron-f one-phonon continuum, which is 
not described by MA^°\ only by MA^^^ and higher lev- 
els [10]. Indeed, as shown in Fig. 2, there is a significant 
difference between the corresponding n — I curves, and 
fnlii'^) does have an onset at wth — ^ even though it be- 
comes hard to see at larger A. The n > 2 curves are much 
less affected, in particular their onset roughly agrees with 
that predicted at MA*^°^ level. Similar behavior is found 
in ID, but the peak in each /n(w) moves towards the 
low-energy threshold [18], as expected since the ID free 
electron density of states is singular at the band-edge. 

In Fig. 3, we plot the first, to our knowledge, com- 
plete set of OC curves reported for a 3D Holstein po- 
laron, using both MA and DMC methods. On the whole, 
the agreement is excellent, especially between MA*^^-' and 
DMC. We find that MA^^^ captures all the qualitative 
features of the full OC, as well as being able to resolve 
finer structure near the absorption onset. In particu- 
lar, the "shoulder" that develops on the low-energy end 
of the OC spectrum could never be captured with per- 
turbational methods. In the asymptotic A ^ 0,oo lim- 
its the curves are nearly indistinguishable, as expected 
since MA becomes exact in these limits. In the crossover 
regime A ~ 1, the differences between MA*^°^ and MA^^^ 
are largest, as are those between MA and DMC results. 
Nevertheless, MA does a good job overall, especially con- 
sidering that it is an efficient analytical approximation. 

The shapes of the OC curves can now be understood 
using Eq. (2) and the data shown in Figs. 1 and 2. For 
small A, Pn=i is dominant and the OC is basically pro- 
portional to /„^i(cj). Indeed, its rich structure is clearly 
visible in the OC, as is its threshold Wth = ^- The n> 2 
terms serve only to alter the high-energy tail. In the small 
polaron limit, however, P„=i — >■ and the OC is domi- 
nated by large n contributions. Since these /n(w) curves 
have similar shapes and are shifted by with respect to 
one other, the OC mirrors the LF Poissonian distribution 
of P„. The peak location is also in good agreement with 
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FIG. 2. (color online) MA*^' functions fi^'{ij), for (a) small 
(b) medium, and (c) large el-ph coupling, at ^l/t = 0.5, = 
0.005 in 3D. Also shown is fi°2i{ui) (dashed line). 



2|i?o| = 12At, expected as A — >■ oo [i^]. This is because 
when the electron is moved to a neighboring site by op- 
tical absorption, it loses the polaronic binding energy Eq 
and it also leaves behind excited phonons with the same 
energy. The structure of our expression (2) proves that 
not only the peak energy, but the very shape of the OC 
curve is determined by the probabilities Pn"^ . In ID (not 
shown), the /n(w) functions are more peaked and indi- 
vidual contributions can be seen in the OC, as f2-spaced 
kinks. The agreement with available ID numerical data 
is of similar quality with that of Fig. 3 [In]. 

Consider now the DMFT-like no- vertex correction ap- 
proximation, in which the two-particle Green's function 
in Eq. (3) is replaced by a convolution of polaron spec- 
tral weights (for details, see the supplementary material). 
The scaled result for A = 1, kgT = O.OIJI, is shown 
in Fig. 3(c). A prominent feature is the peak (marked 
by arrow) below ojth, in agreement with low-T data in 
Ref. ['.)], where it was identified as an excitation from 
the GS into the first bound-state. Indeed, the peak is at 
^ 0.12t, their energy difference. The peak likely vanishes 
at T = 0,7; = 0, however this shows that a naive con- 
sideration of the convolution may lead to a qualitatively 
wrong idea of the OC structure. DMC and MA results 
do not show this peak, although the second bound state 
is visible in their spectral weight (see the supplementary 
material) . Its absence implies a vanishing matrix element 
in Eq. (4), likely due to the different symmetry of the 
polaron wavefunction in the two states [19]. Clearly, Eq. 
(2) properly accounts for such selection rules. More anal- 
ysis is provided in the supplementary material found at 
the end of this article. 

To further prove the accuracy of MA, we compare in 
Table I the total integrated OC, S{oo) = dLU(j{uj), 
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FIG. 3. (color online) OC in units of ctq = neh'^ for a 3D Holstein polaron, calculated using MA^"', MA'^' and DMC. The 
inset in panel (c) shows the absorption onset, and the no-vertex correction OC at kaT = O.OIJI, scaled by 10'^ (thin blue line). 



for MA and DMC. The agreement is again very good, 
particularly at small and large A, where we typically find 
less than 5% difference. At intermediate couplings MA 
accounts for over 85% of the total integrated OC, which 
is quite remarkable for an analytical approximation. 

Knowledge of the total integrated OC can also be used 
to verify that MA satisfies the /-sum rule, given by [20]: 



1 2 

■ e^a_Ekin — —T^ H — S{oo) 

TT TT 



(7) 



where i?kin is the GS polaron kinetic energy and D = 
Tre^a/m* is the Drude weight [2].], m* being the polaron 
effective mass. Both these quantities can be calculated 
from the GS properties of the polaron, using known MA 
methods [ ]. fn all cases reported here, we have verified 
that the total integrated OC matches the expected value 
from the /-sum rule to at least three decimal places. 

In conclusion, we have obtained the first accurate an- 
alytical non-perturbational expression for the T = op- 
tical conductivity of a Holstein polaron. It explains the 
shape of the OC curve by explicitly relating it to the 
statistics of the polaron cloud. We find no absorption 
below the threshold at ^l even if a second bound state 
exists, proving that Eq. (2) accounts for selection rules. 
Moreover, this work can be generalized for OC studies of 
polaronic systems with momentum-dependent el-ph cou- 
pling, and also disorder, in any dimension. 
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TABLE I. Total integrated OC S{oo)/crf, for the parameters 
used in Fig. 3, for MA'''^ MA^^^, and DMC. 
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SUPPLEMENTARY MATERIAL 
Calculation of the matrix elements within MA'"' 

First, we need to calculate the generalized single par- 
ticle Green's functions: 



is the momentum average of the bare propagator. 

On the other hand, from its Lehmann representation: 



K=o\i^a)i'^a\c, 



q, n) 



Ea + ir] 



(14) 



SO its residue for uj = Eq^ the GS energy, is the product 

w) = (-q,n|cqG(w)c|^| - Q,to) (8) (0|Ck=olV'a)(V'al4l " '^)- The first term is related to 

the GS quasiparticle weight Zq, and the second is the 
overlap that we need. It follows that: 



q,n) = -^^e-'i-^'(6l)"|0) (9) 



where 



We use Dyson's identity G{uj) = Go{uj) + G{ijj)VGQ{uj) 
where V is the el-ph interaction and Gq{lj) is the re- 
solvent for Uf) ^ n - V, plus the MA^'') one-site 
cloud restriction, to find i^„„i(q, Q, ti^) = Go(Q,w — 
"iri)[(5q,Q^„,m"-! + "i5/„^„j_i(q, w) + g/„,„i+i(q, w)]. 
Here, the free propagator is Go(k,w) — {uj + irj — 
ek)-\ and /„,„i(q,w) = ]^ I]q ^'nm(q, Q, w) are par- 
tial momentum averages related to the ones calculated in 
Ref. [j ]. They can therefore be calculated similarly; how- 
ever, they depend on q only through Sq, therefore they 
are even functions whose contribution to (T+(a;) vanishes 
after the sum over q because of the odd prefactor sing 
coming from the current operator. As a result, within 
MA*^°^ we can replace 

in Eq. (5). The delta functions remove the sums over 
m, Q. All that is left, then, is to find the residues 
|(V'o|Cq| — q,n)p. To achieve this, consider the gener- 
alized single-electron Green's functions: 



F„(q,a;) = (0|ck=oG'(w)cl,| -q,n). 



(10) 



These are also similar to the single particle Green's func- 
tions calculated in Rcf. [ ], and can be calculated in 
terms of continued fractions by the same means. Within 
MA'^^^ they are equal to: 

f-„(q, c) ^ ~ ) (c.)G(k = 0,co). (11) 

go{uj - nil) 

where G(k, w) = (0|ci^G(a;)c|^|0) is the usual single- 
particle polaron Green's function, and we define 
^ra,fc('^) — An{Lu)An^i{uj) ■ ■ ■ Ak{i-j) whcrc the continu- 
ous fractions are [i, 2]: 



nggoiuj - nfl) 



1 - 550(1^ ~ nD,)An+i{uj) 



and 



£'o(w) = ^^Go(q,w) 



(12) 



(13) 



(V'oi4i 



q, n) 



Go(q, cj — nft) 
ga{uj - nn) 



%A„,i(£;o), (15) 



The ^-independent part can be shown to be related to 
the probability to have n-particles at the electron site in 
the GS. In Rcf. [ ], these probabilities were found to be: 



^\A^AEo)\' 



(16) 



The oj-dependent terms can be grouped together with the 
one from i^„m(q, Q, w) into the corresponding function 
fn\uj), whose expression is given in Eq. (6). 

The MA^^) calculations proceed along similar lines, but 
are much more involved because of the enlarged basis. 
The full details will be presented in a longer publication. 



The approximation of no current vertex corrections 

The essence of this approximation is to replace the two- 
particle Green's functions appearing in cr'^(cj), namely 
('0o|cqCQG(aj)cJjCq|-0o): by a convolution of the single- 
particle spectral weights. The procedure is reviewed in 
detail in Ref. [4] . The regular part of the conductivity, 
at finite T and on a cubic lattice, is then found to be 
[their Eq. (19)]: 

A poo 

a{uj) = -^Y^ sin^ k J duj'Aik, cj')^(k, uj + uj') 

x[f{Lo')-f{oj + uj')]. (17) 

Here A(k,aj) is the polaron spectral weight and /(e) = 
{e^'^ + l)~^ is the Dirac function. This is the same formula 
used in DMFT, see for example Eq. (3) in Ref. [■>]. The 
difference is that in DMFT, the integral over the Brillouin 
zone is replaced by an integral over energies e = ek (they 
also replace uj' — !■ v). This explains the appearance of 
the density of states of the free lattice, and of the 
DMFT vertex which is the equivalent of the sin^ k in 
the above equation. Of course, within DMFT one uses 
the DMFT self-energy in the spectral weight. 

In order to limit the quantitative differences of the 
DOS used in DMFT, which is for a Bethe lattice, whereas 
our other results are for a cubic lattice, we use Eq. (17) 
for a cubic lattice to test this approximation. Also, for 



6 



convenience, we use the MA self-energy in the spectral 
weight. This is simpler because in order to calculate the 
DMFT self-energy, one needs to go through iterations to 
reach self-consistency. In any event, the MA and DMFT 
self-energies are qualitatively and even quantitatively in 
fairly good agreement, as shown in the Appendix of Ref. 
[ ]. Moreover, it has been shown that the MA spectral 



weights are very accurate, satisfying eight sum rules ex- 
actly (within MA^^)), therefore using the MA self-energy 
and its corresponding spectral weight can lead, at most, 
to modest quantitative differences. 

Since the MA*-^) self-energy is momentum-independent 
(momentum dependence appears only from the MA^^^ 
level), we can carry the integral over the Brillouin zone 
as follows. We rewrite the spectral weights as 



111 1 

--G(k,w) = — — — -. — — - 

TT Ztt \_uj + iij — eic — ^{^) UJ — irj — — 1j* [uj) 

The product of two spectral weights will result in a sum of 4 terms that can be further factorized, for example: 



oj' + — ek — S(w') cj + cj' + it] — Ek ^ Ti{uj + uj') 

1 



UJ -T.{uj + uj') + T.{uj') 



uj' + ij] — — ^("^0 UJ + uj' + irj — ei^ — + uj') 



(18) 



As a result, the integral over the Brillouin zone now be- 
comes a simple momentum average of the bare propa- 
gator (with shifted frequency) times the sin^ kx term. 
One integral can be done analytically, and the other two 
numerically like other similar momentum averages cal- 
culated for MA(^\ see Ref. [ ]. One is then left with 
evaluating the integral over uj' numerically as well. 

One can approach the T = limit either directly by 
replacing the Fermi-Dirac distribution with a Heaviside 
function, as done in Ref. [4] ; or by calculating the T — > 
limit of the ratio detailed in Refs. [5, 6]. Either approach 
has its own numerical challenges, and overall we find both 
procedures to be more time consuming than evaluating 
our OC formula, even when the computationally trivial 
MA self-energy is used in the spectral weight. Note, also, 
that our OC formula is formulated directly for T 
and therefore one needs not worry about how to properly 
approach this limit. 

The k = spectral weight for the polaron, at interme- 
diary coupling A = 1, is shown in Fig. 4. The GS peak 
still has considerable weight Zq ^ 0.5, and just above 
it we see the peak corresponding to the second bound 
state, followed by the continuum and higher energy fea- 
tures. The rough expectation is that a convolution of 
this curve with itself will show a first peak when the two 
curves are shifted by Ei — Eq, the energy difference be- 
tween the two bound states; also, absorption into the 
continuum will appear for frequencies uj > fl. 

Thus, the naive expectation based on the convolution 
is that as soon as the second bound state is formed, there 
should be absorption into it, at energies below the thresh- 
old fl. This is also what is shown by the full DMFT calcu- 
lations at very low temperatures, see curve of Fig. 2 and 



its interpretation from Fig. 3 in Ref. [ ], in qualitative 
agreement with our own calculation for ksT = O.OIQ, 
shown in Fig. 3(c) of the main article. 

In reality, very careful analysis suggests that the weight 
of this peak decreases monotonically with T, and there- 
fore at T = it is very likely that this approach also pre- 
dicts no absorption below the threshold at fi, in agree- 
ment with the DMC and MA approaches. Our under- 
standing of the disappearance of this peak is that it comes 
from the fact that the k = convolution is multiplied by 
sin^ k, so strictly speaking there is no contribution to 
the OC from the GS momentum, unless some "thermal 
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FIG. 4. Spectral weight A(k = 0,uj) for the 3D Holstein 
polaron, for A = 1,^ = 0.5t,rj = 10"^. The GS energy is 
Eo — —6.65t, and the first bound state is at Ei = —6.53t. 
The continuum onset is visible at Eo + i}. 
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broadening" is allowed at non-zero temperatures. 

Since the exact DMC calculation finds no such absorp- 
tion into the second bound state at any A, the matrix 
element for this transition must be zero, likely because 
of the different symmetry of the phonon clouds in the 
two states, as discussed in detail for ID systems in Ref. 
[7]. Clearly, our MA-based OC formula captures this 
selection rule properly - note that the generalized single- 
particle Green's functions that appear in it are different 
from the Green's function that gives the spectral weight. 

On the other hand, there is nothing in the spectral 
weight that encodes the symmetry of the wavefunction 
associated with a bound state, for example one cannot 
distinguish the parity of a state by looking at the plot of 
the spectral weight. Thus, it seems to us that generically, 
one should expect the convolution approximation to fail 
to properly account for selection rules. 

In fact, this statement can be tested in ID, where as 
shown in Ref. [7], one may expect higher energy bound 
states, some of which are optically active. In other mod- 
els, it has been suggested that there may even be bound- 
states below the continuum which are not visible in the 
spectral weight, because their quasiparticle weight van- 
ishes for symmetry reasons [~ ] - if such states are opti- 
cally active, absorption associated with them would cer- 
tainly not be predicted by the convolution approxima- 
tion. Such a study will be undertaken elsewhere, and 
should clarify whether the agreement in 3D regarding 
the lack of sub-threshold absorption in the no-current 
vertex approximation is accidental, or it actually has a 
more robust explanation. In any event, this is probably 
somewhat of an academic discussion, since such multiple 
bound states appear only at large effective couplings A 
where the bulk of the optical absorption is at very much 



higher energies, and in measurements the very small sig- 
nal at the threshold may be buried in noise. 

Be that as it may, to us it appears that besides be- 
ing numerically more difhcult to evaluate and less tested 
in terms of its accuracy for sum rules etc, the no-current 
vertex approximation also fails to provide the insight into 
the overall evolution of the OC curves that is afforded 
by our new formula. Moreover, MA has already been 
demonstrated to generalize very successfully to models 
where the coupling depends on the phonon and elec- 
tron momentum, and we expect similar success in mod- 
eling their OC. Such models have strongly momentum- 
dependent self-energies, and are unlikely to be described 
accurately with a local-approximation like DMFT. 
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